Quasiparticle tunneling in the Moore-Read fractional quantum Hall State 
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In fractional quantum Hall systems, quasiparticles of fractional charge can tunnel between the edges at a 
quantum point contact. Such tunneling (or backscattering) processes contribute to charge transport, and pro- 
vide information on both the charge and statistics of the quasiparticles involved. Here we study quasiparticle 
tunneling in the Moore-Read state, in which quasiparticles of charge e/4 (non-Abelian) and e/2 (Abelian) may 
co-exist and both contribute to edge transport. On a disk geometry, we calculate the matrix elements for e/2 
and e/4 quasiholes to tunnel through the bulk of the Moore-Read state, in an attempt to understand their relative 
importance. We find the tunneling amplitude for charge e/2 quasihole is exponentially smaller than that for 
charge e/4 quasihole, and the ratio between them can be (partially) attributed to their charge difference. We 
find that including long-range Coulomb interaction only has a weak effect on the ratio. We discuss briefly the 
relevance of these results to recent tunneling and interferometry experiments at filling factor v = 5/2. 

PACS numbers: 



I. INTRODUCTION 

The fractional quantum Hall effect (FQHE) at filling fac- 
tor v = 5/ 2 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7,8 ' 9 ' 10 has attracted strong interest, due 
to the possibility that it may support non-Abelian quasipar- 
ticles, and their potential application in topological quantum 
computation . 11 ' 12 ' 13 ' 14 Numerical studie s 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 in- 
dicate that the Moore-Read (MR) state 23 or its particle-hole 
conjugate state 24,25 are the most likely candidates to describe 
the v = 5/2 quantum Hall liquid. They both support non- 
Abelian quasiparticle excitations with fractional charge e/4, 
in addition to the Abelian quasiparticle excitation with frac- 
tional charge e/2 of the Laughlin typei 2 ^ 6 . 

Edge excitations in the FQHE can be described at low en- 
ergies by a chiral Luttinger liquid model 27 , and quasiparticle 
tunneling through barriers or constrictions was originally con- 
sidered 28,29 in the case of the Laughlin state. Recently the 
transport properties of the v = 5/2 state through a point con- 
tact have also been considered by a number of authors . 30 ' 31 ' 32 
Experimentally the quasiparticle charge of e/4 has been mea- 
sured in the shot noise 33 and temperature dependence of tun- 
neling conductance^ 4 - The latter also probes the tunneling ex- 
ponent, which is related to the Abelian or non-Abelian nature 
of the state, although a direct probe on the statistics based on 
quasiparticle interference is desired. 

The two point contact Fabry-Perot interferometer was first 
proposed for probing the Abelian statistics 35 and later con- 
sidered for the non-Abelian statistics, 36 ' 37 ' 38 ' 39 ' 40 ' 41,42 ' 43 ' 44,45 In 
this kind of setup, quasiparticles propagating along the edges 
of the sample can tunnel from one edge to the other at the 
constrictions formed in a gated Hall bar. Such tunneling pro- 
cesses lead to interference of the edge current between two 
different tunneling trajectories. It has been used in both in- 
tegeri 6 ^ 4 - 7 - and fractional quantum Hall regimes in the lowest 
Landau leveli 4 ^ 9 - Recently, Willett et al^^ implemented 



such a setup in the first excited Landau level and attempted 
to probe the non-Abelian statistics of the quasiparticles in the 
case of v — 5/2 from the interference pattern. 

The interference pattern at v = 5/2 state is predicted to 
exhibit an even-odd variation 38,39 depending on the parity of 
the number of e /4 quasiparticles in the bulk. This would be 
a direct indication of their non-Abelian nature. In their exper- 
iments, Willett et aL— observed oscillations of the longitu- 
dinal resistance while varying the side gate voltage in their 
interferometer. At low temperatures they observed appar- 
ent Aharanov-Bohm oscillation periods corresponding to e/4 
quasiparticle tunneling for certain gate voltages, and periods 
corresponding to e /2 quasiparticle at other gate voltages. This 
alternation was argued to be due to the non-Abelian nature of 
the e/4 quasiparticles^i^ 2 - consistent with earlier theoretical 
prediction;^ 9 - At higher temperatures e/4 periods disappear 
while e/2 periods persist^ - 

There are two possible origins for the e/2 period in the in- 
terference picture: It may come from the interference of e/2 
quasiparticles, or the intereference of e/4 quasiparticles that 
traverse two laps around the interferometer. It is natural to ex- 
pect that the tunneling of the e/4 quasiparticles is much eas- 
ier than that of e/2 quasiparticles. Therefore, the tunneling 
amplitudes of e /4 quasiparticles should be larger than that of 
the e/2 quasiparticles. On the other hand, e/2 quasiparticles, 
being Abelian (or Laughlin type), involve the charge sector 
only and have much longer coherence length than that of e/4 
quasiparticles. 18 In fact it was predicted 18 that the e/2 inter- 
ference pattern will dominate once the temperature dependent 
coherence length for e/4 quaisparticles becomes shorter than 
the distance between the two point contacts, in agreement with 
more recent experiment 3 

In the present paper, we attempt to shed light on the rel- 
ative importance of e/4 and e/2 quasiparticle tunneling in 
transport experiments involving point contacts. By numeri- 
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cally diagonalizing a special Hamiltonian with three-body in- 
teraction that makes the Moore-Read state the exact ground 
state at half-filling, we explicitly calculate the amplitudes of 
e/4 and e/2 quasiparticles tunneling from one edge to another 
through the Moore-Read bulk state in disk and annulus ge- 
ometries. We find the (bare) tunneling amplitude for charge 
e/2 quasiparticles is exponentially smaller than that for charge 
e/4 quasiparticles, and their ratio can be partially (but not 
completely) attributed to the charge difference. These results 
would allow for a quantitative interpretation of the quasiparti- 
cle interference pattern observed by Willett et alM 

The remainder of the paper is organized as follows. In 
Sec. IU we describe the microscopic model of the 5/2-filling 
fractional quantum Hall liquid on a disk and its ground state 
wave functions with and without a charge e/4 or e/2 quasi- 
hole in the center. We introduce the tunneling potential for the 
quasiholes and outline the scheme of our calculation. We then 
present our main results for the case of short-range interaction 
in Sec. lHI Al in which we compare the different tunneling am- 
plitudes of the charge e/4 and e/2 quasiholes. We also map 
the results from the disk geometry to an experimentally more 
relevant annulus geometry and attempt to obtain the leading 
dependence of the tunneling amplitudes on system size and 
inter-edge distance. We discuss the influence of long-range 
Coulomb interaction in Sec. lIIIBl In Sec.|lV] we summarize 
our results and discuss their relevance to recent interference 
measurement in the 5/2 fractional quantum Hall system. 



II. MODEL AND METHOD 

We start by considering a disk on which a v = 1/2 Moore- 
Read fractional quantum Hall liquid resides. The disk geom- 
etry can support both charge e/4 and e/2 excitations at the 
center, providing us with an opportunity to study their tunnel- 
ing to the edge (see Fig. [TJ. Later in the paper, we will also 
map the geometry to an annulus or a ribbon of electrons, thus 
allowing a closer comparison with realistic experimental situ- 
ations, e.g. in the vicinity of a quantum point contact. Our 
system resembles a multiply connected torus with a strong 
barrier studied prieviously 53 in the context of Laughlin quasi- 
particle tunneling. In the half-filling case, we need to consider 
both Laughlin-type Abelian quasiparticles with characteristic 
charge e/2 and Moore-Read-type non- Abelian quasiparticles 
with charge e/4. 

To study the Moore-Read ground states with and without 
an e/4 or e/2 quasihole at a half filling, we start from a three- 
body interaction H$b 

H 3B = - E SiAVl^Sin-r^Siri-rk)}, (1) 

i<j<k 
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FIG. 1: (Color online) Theoretical setup for a disk of a fractional 
quantum Hall liquid, allowing quasiholes to tunnel through the bulk 
from the center to the edge. 



in the lowest Landau level (LLL) representation, 

*mr(^i, z 2 , z N ) = 
1 



Pf 
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is the exact zero-energy ground state of H^b with the smallest 
total angular momentum Mo = N(2N - 3) /2. In Eq. ©, the 
Pfaffian is defined by 



N/2 
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for an N x N antisymmetric matrix with elements . 

The three-body interaction also generates a series of zero- 
energy states with higher total angular momentum, related 
to edge excitations and bulk quasihole excitations. The N- 
electron Moore-Read ground state with an additional charge 
e/4 quasihole at the origin (so the edge also expand corre- 
spondingly due to a fixed number of electrons) has a wave- 
function 
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This state is a zero-energy state with total angular momentum 
Mq+N/2 in the lowest 2^—1 orbitals (one more than needed 
for the Moore-Read state), but not the only one. To generate 
the unique charge e/4 state, we need to introduce a strong 
repulsive interaction for electrons occupying the lowest two 
orbitals 



Ai? e / 4 = Ac^CiCq" 



CO, 



A 



(5) 



On the other hand, the Moore-Read ground state with a e/2 
quasihole (i.e., a Laughlin quasihole, equivalent to two e/4 
quasiholes fused in the identity channel) at the origin, 



where S is a symmetrizer: Si 23 [/i23] = A23 + /231 + fdi2- 
The A^-electron Pfaffian state proposed by Moore and Read 23 



*Mr( Z 1: Z 2,-.-,ZAt) 



JJzi *mr(^i,^2,— ,zn), (6) 



3 



is the unique zero-energy ground state with total angular mo- 
mentum Mo + N in the lowest 2N — 1 orbitals. 

The Moore-Read state [Eq. ©], together with its quasihole 
states [Eqs. © and ©], can therefore be generated by numer- 
ically diagonalizing the three-body Hamiltonian [Eq. (Q]) with 
Eq. ©, the artificial repulsion to generate an e/4 quasihole] 
in corresponding finite number of orbitals using the Lanczos 
algorithm. The wavefunctions can then be supplied to cal- 
culate the tunneling amplitudes of the quasiholes. The same 
numerical procedure can be used to study the tunneling am- 
plitudes for the more realistic situation with a long-range in- 
teraction, in which case the variational wavefunctions are no 
longer eigenstates of the realistic Hamiltonian. For clarity and 
convenience, we will delay the discussion on how to generate 
realistic ground state and quasihole states in the presence of 
long-range interaction (and their comparison with the varia- 
tional states) to Sec. lIIIBI 

To study the tunneling amplitudes of the quasiholes, let us 
first consider a single-particle picture, which will help us un- 
derstand our approach and, later, our results as well. In the 
disk geometry, the single-particle eigenstates are 



\m) = <j> m {z) = {2n2 m m\) 
We assume a single-particle tunneling potential 

V tunneX (9) = V t 5(9), 
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(7) 



(8) 



which breaks the rotational symmetry. Here we calculate the 
matrix element of (fc|Vtunnel(0)|O> related to the tunneling of 
an electron from state \l) to state \k). One can visualize the 
tunneling process as a path along the polar angle 8 = be- 
tween the two states centered around their maximum ampli- 
tude st\z\ = y2l and | jar j = \/2k, respectively. One readily 
obtains 



«p(M)= WunnclWIO 



V t r (jg + 1) 

2tt JW\ 



(9) 



The interesting limit is that we let k and I tend to infinity, 
but keep the tunneling distance fixed at d, i.e., \k — l\ ~ 
\/2k(d/lB) <C (k + I). Alternatively, we can understand d 
through the angular momentum change l 2 B \k — l\/R, where 
R ~ \/2kl b is the azimuthal size of the single-particle state 
with momentum k (or I in this limit). We can show (see Ap- 
pendixfA}, in this limit, 

Vp{ki |) „ YL e -(k-rf/Mk + i+2) „ ^ e -dV(2i B) ^ (1Q) 
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which reflects the overlap of the two Gaussians separated by 
a distance d. 

For quasiparticle tunneling at filling fraction v = 5/2, one 
should, in principle, use wavefunctions in the first excited 
Landau level (ILL). Evaluating the tunneling matrix element 
in the ILL, we obtain an additional prefactor, so 



vf L {k,l) = 



1 



(k-l) 2 
2(k + l) 



v p (k, I). 



(11) 



The sign change in the prefactor at d ~ Is can, unfortunately, 
cause severe finite-size effect for the numerically accessible 
range. Nevertheless, in the thermodynamic limit, the prefactor 
can be approximated by — (fc — l) 2 /2(k + I) ~ — d 2 /(2l 2 B ) 
and, therefore, the leading decaying behavior is essentially the 
same. So we will continue to work in the LLL but expect that 
the leading scaling behavior is the same as in the ILL. 

In the many-body case, we write the tunneling operator as 
the sum of the single-particle operators, 



tunnel V 



We are now ready to calculate the tunneling amplitudes 



T e / 4 = (*mr|T|*£&) and r e / 2 



*mr|T|* 



e/2, 
MR/ 



for 

e/4 and e/2 quasiholes, respectively. For convenience, we 
will set Vt = 1 as the unit of the tunneling amplitudes in 
the following text and figures. As explained in Ref. HH 
the matrix elements consist of contributions from the respec- 
tive Slater-determinant components \li, Z/v) £ vPmr and 

\ki, fcjv) £ ^mr or ^mr- Non-zero contributions en- 
ters only when \l±, In) and \k\, kpj) are identical except 
for a single pair / and k with angular momentum difference 
k — I = N/2 or N for the quasihole with charge e/4 or e/2. 

For clarity, we also include a pedagogical illustration of the 
procedure for calculating the tunneling matrix elements in the 
smallest possible system of four electrons in Appendix IE1 



HI. RESULTS 
A. Short-range interaction 

Systems of up to six electrons can be worked out pedagog- 
ically using Mathematica as illustrated in Appendix [B] For 
larger systems, we obtain the exact Moore-Read and quasihole 
wavefunctions by the exact diagonalization of the three-body 
Hamiltonian [Eq. (fl~|i] using the Lanczos algorithm. The tun- 
neling amplitudes are then evaluated as explained in Sec. UT1 
Figure[2fa) plots the tunneling amplitudes for the e/4 and e/2 
quasiholes in the Moore-Read state as a function of electron 
number. The result for the e/4 quasihole shows a weak in- 
crease for N < 10 followed by a decrease for N > 10. On 
the other hand, the result for the e/2 quasihole shows a mono- 
tonic decrease as the number of electrons increases up to 14. 
In the largest system, the ratio of the two tunneling matrix el- 
ements is slightly less than 20. For comparison, we also plot 
the tunneling amplitudes T e / 3 and T 2e / 3 for the e/3 and 2e/3 
quasiholes in a Laughlin state at v = 1/3 in Fig.|2jb). We also 
observe a bump in the tunneling amplitude T 6 / 3 for charge 
e/3, followed by a monotonic decrease. We thus expect the 
r e / 4 would eventually also show a monotonic decrease for 
large enough systems. T 2e / 3 for charge 2e/3 shows a much 
faster decrease, consistent with its larger charge, and thus a 
larger momentum transfer for the same tunneling distance. 

Due to the finite-size bumps in T e / 4 and T e / 3 for charge 
e/4 and e/3, it is difficult to extract the asymptotic behavior 
in the tunneling amplitudes for these quasiparticles. However, 
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FIG. 2: (Color online) tunneling amplitude as a function of number 
of electrons for (a) e/4 and e/2 quasiholes in the Moore-Read state 
at half filling and (b) e/3 and 2e/3 quasiholes in the Laughlin state 
at is = 1/3. 
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FIG. 3: (Color online) The ratio of tunneling matrix elements for e/2 
quasiholes to e/4 quasiholes in the Moore-Read state at half filling, 
and for 2e/3 quasiholes to e/3 quasiholes in the Laughlin state at 
1/3 filling as a function of number of electrons. 



we may expect that such finite-size corrections also exist in the 
tunneling amplitude for charge e/2 and 2e/3 so we can extract 
the asymptotic behavior in their ratios. Fortunately, this is 
indeed the case. We plot r e / 2 /r e/4 and r 2e / 3 /r e / 3 in Fig.H 
We find the ratios can be fitted very well by exponentially 
decaying trends for almost all finite system sizes. The fitting 
results are 



T e/2 /T e/4 „ L78e -0.25W 
r 2e/3/ r e/3 _ q ^ 



-0.40^ 



(13) 
(14) 



As will be discussed later, the exponents are related to the 
charge of the quasiholes and, to a lesser extent, to corrections 
due to sample geometry, perhaps also to the influence of the 
neutral component of the charge-e/4 quasiparticles. Qualita- 
tively, the constant in the exponent of the ratio r e / 2 /r e / 4 is 
found to be smaller than that for r 2e / 3 /r e / 3 , consistent with 
the smaller charge and thus smaller charge difference in the 
half-filled case. 

One may question that the tunneling amplitude for a quasi- 
hole from the disk center to the disk edge may be different 
from that for edge to edge, as in the realistic experimental 



FIG. 4: (Color online) (a) The tunneling amplitude T e ^ 4 and (b) 
the tunneling amplitudes r e//2 as functions of edge-to-edge distance 
d(n,N). Data is shown up to n = 100 quasiholes and N = 14 
electrons. The data points at d = or n — > oo limit are exact results 
as explained in Appendix [Cl 



situations. In particular, the former can contain a geometric 
factor, which can be corrected by mapping the disk to a an- 
nulus (or a ribbon) by inserting a large number of quasiholes 
at the disk center, from which electrons are repelled (see Ap- 
pendix ICl for technical details). Inserting n quasiholes to the 
center of a disk of N electrons in the Moore-Read state, we 
can write the new wavefunction as 



ne/2 
MR 



' N 

n 



MR, 



(15) 



so that each component Slater determinant gets shifted into a 
new one to be normalized. The first n orbitals from the center 
are now completely empty and the electrons are occupying 
orbitals from n to n + 2N — 3. This transformation, of course, 
also changes the tunneling distance to 



d{n, N) jl B = y/2{n + 2N -2) - y/2n. 



(16) 



So we can plot data using d(n, N), rather than n. Similarly, 
we can make the same transformation for the Moore-Read 
state with either an additional charge e /4 excitation or an ad- 
ditional charge e /2 excitation at the inner edge defined by the 
inserted n quasiholes. Thus, we can calculate the tunneling 
amplitudes under the mapping from disk to annulus. 

In Fig.lH we show the tunneling amplitudes T e / 4 and T e / 2 
for up to n = 100 quasiholes. We plot them as functions 
of tunneling distance d, which decreases as n increases. It 
is interesting to note that finite-size effects diminish beyond 
d > 61b for charge e/4 and d > 51b for charge e/2. For 
comparison, we plot the ratio rW 2 /r e / 4 as a function of d in 
Fig- El We find that, when we insert more than one quasihole, 
the ratio of the tunneling amplitudes falls onto a single curve, 
regardless of the system size N and the number of quasiholes 
n. The curve can be fit roughly to 



r e / 2 (d) 



-0.083(d/Z B ) 2 



(17) 



r e / 4 (d) 

We point out that a few points in Fig. [5] can been seen de- 
viated from this behavior. They correspond to the largest d 
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FIG. 5: (Color online) The ratio of tunneling amplitudes r e/2 /r e/4 
as a function of edge-to-edge distance d(n, N). We also plot Eq. l !17t 
as the solid line to guide the eye. The dashed line is the theoretical 
estimate based on the charge component only [Eq. (19)]. 

for a given N, meaning that there is no quasihole in the bulk, 
thus corrspond to the bulk-to-edge instead of the edge-to-edge 
tunneling. 

It is worth pointing out that such behavior is not completely 
unexpected; in fact, it reflects the asymptotic behavior of the 
single-particle tunneling matrix and the corresponding charge 
of the quasiparticles. To see this, we note that for a charge 
q quasihole to tunnel a distance of d, one electron (in each 
Slater determinant) must hop by a distance of qd/e for the ex- 
act momentum transfer. According to the asymptotic behavior 
in Eq. ( fTUl ), we expect 

Therefore, we expect 

_ e -[(rf/2) 2 -(d/4) 2 ]/(2/ B ) 2 = e -0.047(d/i B ) 2 (19) 
pe/4 ' v ' 

which we also include in Fig.|5]for comparison. 

We thus find that both variational wavefunction calcula- 
tion and qualitative analysis suggest that the tunneling am- 
plitude of the e/2 quasiparticles is smaller than that of the e/4 
quasiparticles by a Gaussian factor in edge-to-edge distance 
d, which is the main results of this paper. There is, however, a 
quantitatively discrepancy in the length scale associated with 
the Gaussian dependence between Eqs. ( fT9] l and illi . This 
indicates that the Gaussian factor in single-electron tunnel- 
ing matrix element only partially accounts for the Gaussian 
dependence; the remaining decaying factor thus must be of 
many-body origin, whose nature is not clear at present and 
warrants further study. 

B. Long-range interaction 

So far, we have discussed the tunneling amplitudes us- 
ing the variational wavefunctions, which are exact ground 
states of the three-body Hamiltonian. These wavefunctions 
are unique, but in general not the exact ground states of any 



generic Hamiltonian one may encounter in a realistic sample. 
In reality, long-range Coulomb interaction is overwhelming, 
although Landau level mixing can generate effective three- 
body interaction. 54 In this subsection, we explore the quasi- 
hole tunneling in the presence of long-range Coulomb inter- 
action. The central questions are the following. First, how 
can we generate both non-Abelian and Abelian quasiholes 
in practice? Remember now we do not have the variational 
Moore-Read state as the exact ground state, so the variational 
quasihole states are also less meaningful. We attempt to gen- 
erate and localize quasiholes with a single-body impurity po- 
tential; then how close are the corresponding wavefunctions to 
the variational ones? Second, suppose we have well-defined 
quasihole wavefunctions, are the results on the tunneling am- 
plitude obtained in the short-range three-body interaction case 
robust in the presence of long-range Coulomb interaction? 

For a smooth interpolation between the short- and long- 
range cases, we introduce a mixed Hamitonian 

H x = (l-\)H C + \H 3B , (20) 

as explained in our earlier works. Here, the dimension- 
less A interpolates smoothly between the limiting cases of 
the three-body Hamiltonian H-sb (A = 1) and a two-body 
Coulomb Hamiltonian He (A = 0). H c also includes a back- 
ground confining potential arising from neutralizing back- 
ground charg e di stributed uniformly on a parallel disk of ra- 
dius R — V4JV, located at a distance D above the 2DEG. 
Using the symmetric gauge, we can write down the Hamilto- 
nian for electrons in the ILL as 

H C = \Y1 V in4n+l C i c n+lCm + ^ U ™ C L C ™> ( 21 ) 
mnl m 

where c] n is the electron creation operator for the first ex- 
cited Landau level (ILL) single electron state with angular 
momentum m. V^'s are the corresponding matrix elements 
of Coulomb interaction for the symmetric gauge, and t/ m 's 
the corresponding matrix elements of the confining potential. 
We choose D = 0.6 so the ground state can be well described 
by the Moore-Read state. 

To be experimentally relevant, we also want to generate the 
quasihole states by a generic impurity potential, rather than by 
the artificial interaction [Eq. Q] we used above to generate 
the unique e/4 quasihole state in the three-body case. We 
consider a Gaussian impurity potential, 55 

H imp (W,a) = Wj2^ m2/2s2c l c m, (22) 

m 

which will trap at the disk center an e/4 or e/2 quasihole de- 
pending on its strength^ Here, s characterizes the range of 
the potential. Note Hi mp = WcqCq is the short-range limit 
(s — ► 0) of the Gaussian potential in Eq. (1221 . W is always 
expressed in units of e 2 / (cIb)- 

Earlier studies 18 ' 55 have identified s = 2.0 as a suitable 
width for the Gaussian trapping potential, which is of roughly 
the radial size of a quasihole. So we use this value exclu- 
sively in the following discussion. One expects that for small 
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FIG. 6: (Color online) (a) Energies of the e/4 and e/2 quasihole 
states E^ h , measured from the corresponding ground state in the 
momentum Mq subspace E x , as a function of the strength W of 
the Gaussian trapping potential at the disk center with s — 2.01b for 
the v — 5/2 state with a mixed Hamiltonian H\ (A = 0.5). The e/4 
quasihole state is energetically favorable for 0.032 < W < 0.137. 
(b) Overlaps of the e/4 and e/2 quasihole states l^ 9 ^ 1 ) for the 
mixed Hamiltonian with the corresponding variational states I^mk) 
[Eqs. Hi and ©]. 



W, the system remains in the Moore-Read phase without any 
quasihole excitation in the bulk; for later reference, we use E x 
to denote the ground state energy in the momentum subspace 
of M = M = N(2N - 3)/2. As W increases, the im- 
purity potential first tends to attract a charge-e/4 quasihole, 
the smallest charge excitation, at the disk center. This would 
be reflected in the sudden angular momentum change from 
Mq to Mo + N/2 of the global ground state, which is also 
characterized by a depletion of 1 /4 of an electron in the elec- 
tron occupation number at orbitals with small momentum. We 

e/4 

use E x to denote the ground state energy in the subspace of 
M = Mq + N/2. When W is increased further, one can 
trap a charge-e/2 quasihole at the center, with ground state 
having the total angular momentum of Mo + N; in this mo- 

e/2 

mentum subspace, we use E x to denote the ground state 
energy. We illustrate this scenario for the case for A = 0.5 
in Fig. [6f a), in which we plot the energies of the e/4 and 
e/2 quasihole states E e x ^ and E e x 2 , measured from the corre- 
sponding E x . More precisely, the e /4 quasihole state is ener- 
getically favorable for 0.032 < W < 0.137. At W < 0.032, 
we find E° x < E c x /4 < E e x /2 , while at W > 0.137, we find 
El > E e x /A > E e x /2 . 

To understand how good these wavefunctions are, we plot, 
in Fig.|6jb), the overlap of the e/4 quasihole state I'l'^/ 4 ) with 
the corresponding variational state l^iT?) P^l- ©1> as we ^ 
as the overlap of the e/2 quasihole state \^ e x 2 ) with the cor- 
responding variational state l^/fl) [Eq. ©]. We find that for 
intermediate W the overlaps are larger than 97%. At small 
W, l^^ 2 ) does not agree with |4'^{ 2 ^) well, but they have 
very large overlap when the charge-e /2 quasihole state energy 

E e J 2 is lower than the correponding energy E x of the Moore- 
Read-like state. On the other hand, it is a little surprising 



FIG. 7: (Color online) The ratio of the tunneling amplitudes for a 
mixed Hamiltonian as a function of number of electrons. The mix 
parameter A = 0.5 and the width and strength of Gaussian potential 
are W = 0.1 and s — 2.0, respectively. The background potential is 
located at the distance D — 0.61b- The dotted line is the exponen- 
tially decaying trend line [Eq. J13H for the pure three-body case as 
shown in Fig. [3] 



to see the excellent agreement between \^ x ) and l^ffl)' 
as they are generated by different Hamiltonian [Eqs. (|20T > and 
([TJ] with different trapping potential [Eqs. d22l and Q] re- 
spectively. But we note that the smooth Gaussian trapping 
potential does favor the charge-e/4 quasihole state, which has 
no simultaneous occupation of the lowest two orbitals. 

Therefore, we expect that with a moderate mixture of the 
long-range Coulomb interaction, the results on the tunnel- 
ing amplitudes are rather robust. In particular, we choose 
W = 0.1, at which we have E° x > E e x /2 > E e x /4 and 

at which both |<< 4 fl l*A /4 )| 2 ^d \(*% 2 R \*f)\ 3 are very 
close to 1 . For example, we plot the ratio of tunneling ampli- 
tudes r e / 2 /r e / 4 as a function of the number of electrons for 
A = 0.5 in Fig. [7] The data points are in good agreement with 
the trend [Eq. ( fT3] )l obtained earlier for the pure three-body 
case. 

So far, we have shown a case where the presence of the 
long-range interaction has very weak effects on the results of 
tunneling amplitudes. However, in general, one can expect 
such an agreement becomes worse as one move farther away 
from the pure repulsive three-body interaction in the parame- 
ter space. To present a more quantitative picture, we plot in 
Fig. [8] the ratio of tunneling amplitudes r e / 2 /r e / 4 (without 
inserting quasiholes at the center, i.e. n = 0) as a function of 
A for the 12-electron system with the mixed Hamiltonian H\ 
and a Gaussian trapping potential (W — 2.0, s = 1.0). The 
ratio remains as a constant from A = 1 down to 0.2, before 
it fluctuates significantly; the fluctuation is believed to be re- 
lated to the stripe-like phase near the pure Coulomb case in 
finite systems, as also revealed in earlier work. 18 We point out 
that recent numerical work suggests that the spin-polarized 
Coulomb ground state at v = 5/2 is adiabatically connected 
with the Moore-Read wave function for systems on the surface 
of a sphere, 22 so the large deviation may well be a finite-size 
artifact. 

Varying parameters, such as W, s, and D, can also lead 
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FIG. 8: (Color online) The ratio of tunneling amplitudes r e/2 /r e/4 
as a function of the mixing parameter between the three-body inter- 
action and Coulomb interaction in a 12-electron system at half-filling 
in case of D — QMb- The width and strength of Gaussian potential 
are W = 0.1 and s = 2.0 respectively. The green dot at A = 1 is 
the value for the short-range three-body interaction case obtained in 
Sec. lmAl 



to larger deviation from the pure three-body case, although 
we find in generic cases r e / 2 /r e/4 remains small. We re- 
mind the reader that the Moore-Read phase is extremely frag- 
ile. Therefore, we have rather strong constraints on param- 
eters when both the Moore-Read-like ground state and the 
quasihole states should subsequently be good description of 
the ground states as the impurity potential strength increases. 
For example, the window of D for the ground state at W = 
to be of Moore-Read nature in the pure Coulomb case is very 
narrow (0.51 < D/l B < 0.76 for 12 electrons in 22 orbits^); 
in this range, the effect of the background potential parame- 
ter D on the ratio of tunneling amplitudes is negligible (less 
than 1% variation). Therefore, we conclude that the small ra- 
tio of r e / 2 /r e / 4 is robust in the presence of the long-range 
interaction as long as the system remains in the Moore-Read 
phase. 



IV. DISCUSSION 

In this work we use a simple microscopic model to study 
quasiparticle tunneling between two fractional quantum Hall 
edges. We find the tunneling amplitude ratio of quasiparticles 
with different charges decays with a Gaussian tail as edge-to- 
edge distance increases. The characteristic length scale as- 
sociated with this dependence can be partially accounted for 
by the difference in the charges of the corresponding quasi- 
particles. More specifically, we find the tunneling amplitude 
for a charge e/4 quasiparticle is significantly larger than a 
charge e/2 quasiparticle in the Moore-Read quantum Hall 
state, which may describe the observed fractional quantum 
Hall effect at the filling factor v = 5/2. This result was antic- 
ipated in Ref. HH in which the authors outlined a microscopic 
calculations that is similar to the discussion in Sec. lIII Al (see 
their Appendix B). 

It is worth emphasizing that what we have calculated here 
are the bare tunneling amplitudes. Under renormalization 



group (RG) transformations, both amplitudes will grow as 
one goes to lower energy/temperature, as they are both rel- 
evant couplings in the RG sense. The ratio between them, 
r e / 2 /r e / 4 , will decrease under RG, because T e / 4 is more rel- 
evant than r e / 2 , which renders T e/2 even less important than 
r e / 4 at low temperatures. This is clearly consistent with tun- 
neling experiments involving a single point contact, 33,34 where 
only signatures of e/4 tunneling is seen. 

However, the importance of the two kinds of quasiparti- 
cles can be reversed in interferometry experiments that look 
for signatures from interference between two point contacts. 
This is because the interference signal depends not only on 
the quasiparticle tunneling amplitudes, but also their coher- 
ence lengths when propagating along the edge of fractional 
quantum Hall samples. Recently, Bishara and Nayak 41 found 
that in a double point-contact interferometer, the oscillating 
part of the current for charge q quasiparticles can be written 
as 



z^oc 7 |rW||rW| e -i-K 



cos ' 



where v[ % are the charge q quasiparticle tunneling amplitudes 
at the two quantum point contacts 1 and 2 with a distance of 
xi2- 7 is a suppression factor resulting from the possible non- 
Abelian statistics of the quasiparticles. For q = e/2, we have 
s = 1, while for q = e/4, 7 = ±l/\/2 (or 0) when we 
have even (or odd) number of e/4 quasiparticles in bulk. The 
sign depends on whether the even number of e/4 quasiparti- 
cles fuse into the identity channel (+) or the fermionic chan- 
nel (— ). $ is the flux enclosed in the interference loop and 
$0 = hc/e is the magnetic flux quantum. The phase S^ qS> is 
the statistical phase due to the existence of bulk quasiparticles 
inside the loop and a the phase arg^iT^). At a finite tem- 
perature T, the decoherence length for the quasiparticle 
in the Moore-Read state is 41 



(23) 



r(9) _ 



1 



' (9) 



9 



2ttT 



(24) 



where v Cj7l are the charge and neutral edge mode velocities 
and g'fX the charge and neutral sector scaling exponents for 
charge q quasiparticles, respectively. Earlier studies by the 
authors 17,18 found that the neutral velocity can be significantly 
smaller (by a factor of 10) than the charge velocity, leading to 
a shorter coherence length L^ 4 * 1 for charge e/4 quasiparti- 
cles (less than 1/3 of L^ 2 ^ for charge e/2 quasiparticles in 
the Moore-Read case, as L^ 2 ^ depends on v c only because 
it is Abelian and gn^ 2 ^ = 0). Finite-size numerical analysis 56 
maps out the dependence of L^ 4 ^ and L^ 2 ^ on the strength 
of confining potential parametrized by D for the Moore-Read 
state, as summarized in Fig. [9] for T = 25 mK used in the 
recent experimental study. 50 Depending on the size of the in- 
terference loop and the strength of the confining potential, the 
edge transport may exhibit both e/4 and e/2 quasiparticle in- 
terference, e/2 quasiparticle interference only, or no quasi- 
particle interference. In particular, the observation of the e/4 
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FIG. 9: (Color online) Decoherence length L$ as a function of D for 
both e/4 (upper line) and e/2 (lower line) quasipaticles in the Moore- 
Read Pfaffian state. We choose a temperature T = 25 mK to al- 
low a direct comparison with experiment. 50 The broken lines above 
D — 0.621b are obtained by extrapolation, as the Moore-Read-like 
ground state is no longer stable in a system of 12 electrons in 26 
orbitals. We note a stripe phase may emerge below D = /s.— 



quasiparticle interference depends sensitively on the length of 
the interference loop due to the effect of the confining poten- 
tial strength on the neutral mode velocity v n . 

We close by stating that by combining the small ratio be- 
tween e/4 and e/2 quasiparticle tunneling matrix elements 
and the fact that e/2 quasiparticle has longer coherence length 
along the edge, it is possible to provide a consistent inter- 
pretation of the recent tunnelin g 33 ' 34 and interference 5 ^! ex- 
periments. Similar conclusions were reached in a recent 
comprehensive analysis 52 of the interference experimentsi^i 
We would like to caution though that a complete quantita- 
tive understanding of the experiments is not yet available at 
this stage due to our incomplete understanding of the actual 
ground state, mesoscopic effects, and the possible oversimpli- 
fications of the microscopic model (for example edge recon- 
structio n 57 ' 58 ' 59 may occur and complicate the analysis signif- 
icantly). Nonetheless, we hope that the quantitative analysis 
presented here can help solve the puzzle. 



APPENDIX A: SINGLE-PARTICLE TUNNELING MATRIX 
ELEMENTS IN THE LARGE DISTANCE LIMIT 

In the disk geometry, the single-particle eigenstates are 

|m) = 4> m (z) = (27r2 m m!)- 1 / 2 z m e~ |z|2/4 . (Al) 
If we assume a single-particle tunneling potential 

funnel = V t 5(9), (A2) 

the matrix element of (fc|V t |Z), related to the tunneling of an 
electron from state \l) to state |fc), is 



V p (k,l) = (fc|Vtunncl|0 



2n 



(A3) 



Using beta functions 



B{x,y) 



T(x)T(y) 



(A4) 



T(x + y) ' 

we can rewrite the dimensionless tunneling matrix element as 

1/2 



v P (k, I) 



2lTV p (k, I) 

Vt 



fl (fc±I + 1) fc±I + i) 



(A5) 

We are interested in the limit of large I and large k, where we 
can use the asymptotic formula of Stirling's approximation 



B{x,y) ~ V2tt 



x x-l/2y V -l/2 



i \x+y-l/2 

(x + y) ™ 1 
for large x and large y. Therefore, we have 



(A6) 



v p (k,l) 



(¥ + i) 



(k + i) k+1/2 (i + if +1/2 

For convenience, we define 



(A7) 
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S = and M = . 

2 2 

If we further take the limit of S ^> \M\, we find 



v p (k,l) 



1 



M 
5+1 

M 
5 + 1 



M-=(S + l/2) _ M 2 
g 2(S + l) a g 

_ M 2 
g 2(S + 1) _ 



S + M 1 
2 4 



1 



M 
5 + 1 



M 



(A8) 



5+1 



S — M 1 
2 4 



M 
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APPENDIX B: TUNNELING MATRIX ELEMENTS IN 
FEW-ELECTRON SYSTEMS 

In this appendix, we first illustrate the calculation of the 
tunneling matrix elements in a four-electron system for the 
Moore-Read state. In this case, the normalized Moore-Read 
wavefunction can be written as a sum of Slater determinants 

as 



MR 



10|011110) - >/2| 101101) + |110011) 



13 



(Bl) 



where the ket notation denotes a Slater determinant with 
electrons occuping the single-particle orbitals labeled by 1. 
For example, 1 110011) means the normalized antisymmetric 
wavefunction of four electrons occupying the orbitals with an- 
gular momentum 0,1,5,6 (reading from left to right in the ket). 
This can be obtained by explicitly expanding the Moore-Read 
state with four electrons (trivial with the help of Mathemat- 
ica). The corresponding e/4 quasihole state, similarly, can be 
written as 

^mr = ^=(V3|1010101)-4V6|1001110) 

+ 8V2|0110110) -4V3|0111001)V (B2) 



and the e/2 quasihole state 



e/2 
MR 



1 



(10|0011110) 



3V13 
2V3|0101101 



lOllOOli; 



(B3) 



In the many-body case, we write the tunneling operator as 
the sum of the single-particle operators, 



(B4) 



and 



calculate 

e/4 N 



amplitudes T e / 4 

e/2\ 



the tunneling 

(*mr|T|*^) and r e / 2 = (*mr|T|*^) for e/4 
and e /2 quasiholes, respectively. The matrix elements consist 
of contributions from the respective Slater-determinant 



components Zjv) € $mr and fcjv) € 



or * 



e/2 
MR- 



74 
MR. 



There are non-zero contributions only when the 



two sets {Zi, Z/v} an< ! •■•> fc.zv} are identical except 
for a single pair I and k with angular momentum difference 
k — I = N/2 or N for the quasihole with charge e/4 or e/2. 
One should also pay proper attention to fermionic signs. 
With some algebra, one obtains, for the four-electron case, 

1 



MR 



e/4, 
MR/ 



16V5V3.5) 
5V143 L py 

4^30^(4,6) + 8V3« P (2, 4) 
8V2« p (0,2) + 4%/6v p (l,3) 



and 



MR 



■•/2 ) _ 
MR/ — 



10 



1 

39 

2^30^(2,6) 



(B5) 



Wv p (l,5) + 10V2w p (0,4) 

(B6) 



where, as before, we define 



v p (k,l) 



v t r 



' k+l 
v 2 



27T 



(B7) 



The numerical values for the two tunneling matrix elements 
are 0.213 and 0. 123, respectively, in units of Vt. Therefore, in 
the smallest nontrivial system, we find that the tunneling am- 
plitude for e /4 quasiholes is roughly twice as large as that for 
e/2 quasiholes. The example of the four-electron case illus- 
trates how the tunneling amplitudes can be computed. The 
results are, however, not particularly meaningful as the sys- 
tem size is so small that one cannot really distinguish bulk 
from edge. 

A similar analysis can be performed for a system of six 
electrons with the help of Mathematica. Due to larger Hilbert 
space, we will not explicitly write down the decomposition of 
the ground states and quasihole states by Slater determinants. 
Instead, we only point out that the tunneling matrix elements 
are given by 



MR. 



e/4, 
MR/ 



and 



in units of Vt ■ 



MR 



he/2, 
MR / 



0.267 



= 0.105, 



(B8) 



(B9) 



APPENDIX C: MAPPING FROM DISK TO ANNULUS 

Microscopic quantum Hall calculations are commonly 
based on one of the following geometries (or topologies): 
torus, sphere, annulus (or cylinder), and disk. In a specific cal- 
culation, they are chosen either for convenience, or the need 
for having different numbers of edge(s). On the other hand 
one can also map one geometry to another by means of, e.g., 
quasihole insertion. Here, to connect the theoretical analysis 
with experiment, we perform a mapping from the disk to the 
annulus geometry by inserting a large number of quasiholes 
at the center of the disk, effectively creating an inner edge, as 
the electron density in the center is suppressed by inserting a 
small disk of Laughlin quasihole liquid. 

After inserting n charge e/2 Laughlin quasiholes to the cen- 
ter of an TY-electron Moore-Read state, the ground state can 
be written as 



ne/2 
MR 



' N 

n 



MR, 



(CI) 



where the additional factor transforms each Slater determinant 
into a new one to be normalized. Let us use the case of four 
electrons as in Appendix iBl to illustrate. We note, a Slater de- 
terminant 1 01 11 10) with an addition of n Laughlin quasiholes 
evolves into another Slater determinant |(0™)011110), mean- 
ing that the m-th (in this example, m — 1-4) single -particle 
orbital is now mapped to the (m + n)-th orbital. Due to the 
difference in normalization, the latter determinant should be 
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multiplied by a factor of F(n; 1, 2, 3, 4) with a general form 
of 



where the normalization factor C is, explicitly, 



TV 



F(n; mi ,m 2 ,---,m N ) = 2 nN/2 ]J J^±^. (C2) 



Therefore, when we express Eq. dCU explicitly for Eq. ( IB lb , 
we have 



C 



1 10 



1 



1 1 



1!2!3!4! 13 0!2!3!5! 13 0!1!4!5! 13' 



(C5) 



ne/2 
MR 



FUr. !.2.:!.4i^£|(0")011110) 



/2 

F(n;0,2,3,5)^|(0 n )101101) 
v 13 



+ F(n;0,l,4,5)^|(0")110011) 
v 13 



, (C3) 



where M is a numerical normalization factor. For n = 1 we 
thus obtain exactly Eq. (IB3b as expected. Interestingly, in the 
n — > oo (ring) limit, the normalized wavefunction becomes, 
aymptotically, 



MR — ° 



' =|(0")011110) 



1!2!3!4! VT3 



1 V2 
0!2!3!5!7l3 



|(0™)101101) 



1 1 
0!1!4!5! 



|(0")110011) 



, (C4) 



In the limit of n ^> N, we have v p (n + mi, n + m^) — » 
Vf/(27r). Therefore, the tunneling matrix between the states 
with n quasiholes and n + 1 quasiholes becomes 



<*MRV|*MR 2+e/2 > 



2tt 



1 10 
1!2!3!4! 13 



1 VTo I 1 \/2 



1 



1!2!3!4! y!3 V 0!2!3!5! v!3 

" V2 / i vTo 



0!2!3!5! Vl3 V l!2!3!4! JJz 



(C6) 



The tunneling of e/4 quasiholes can be worked out 
in a similar fashion and we can obtain generically 

T ne/2 | | , ne/2+e/A 
MR, I 1 I * MR 
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